NASA TECHNICAL 
MEMORANDUM 

NASA TM X-53857 


ACCELEROMETER CALIBRATION IN THE LOW g 
RANGE BY MEANS OF MASS ATTRACTION 


By Konrad Reinel 
Astrionics Laboratory 

July 28, 1969 




NASA 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 


MSFC - Form 3190 (September 1968) 


TECHNICAL REPORT STANDARD TITLE PAGE 

2. GOVERNMENT ACCESSION NO. I 3. RECIPIENT'S CATALOG NO. j 


1 . REPORT NO. 1 2. GOVERNMENT ACCESSION NO. I 3. RECIPIENT'S CA 

TMX-53857 


4. TITLE AND SUBTITLE 5 REPORT DATE 

Accelerometer Calibration in the Low g Range ' mm or i q R q 
by Means of Mass Attraction 


6 . PERFORMING ORGANIZATION CODE 


7. AUTHOR(S) 


Konrad Reinel 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


8 . PERFORMING ORGANIZATION REPORT # 


10. WORK UNIT NO. 


11 . CONTRACT OR GRANT NO. 


12. SPONSORING AGENCY NAME AND ADDRESS 


|13. TYPE OF REPORT & PERIOD COVERED 


Technical Memorandum 


14. SPONSORING AGENCY CODE 


15. SUPPLEMENTARY NOTES 


Prepared by Astrionics Laboratory, Science and Engineering Directorate 


16. ABSTRACT 

Mass attraction is used as an equivalent acceleration input to calibrate an accelerometer. 
The upper limit of the acceleration by a reasonable mass size is 10 -9 g in orbit and 10~ 7 g in the 
laboratory. The calibration has been carried out in the laboratory for an electrostatic suspended 
single-axis accelerometer (MESA) with a variable mass attraction. The mass attraction input 
to the accelerometer was a sine wave with the amplitude of 23 nano g. The response of the 
accelerometer to this acceleration input by mass attraction was obtained by data analysis and 
data reduction with a computer. The results of the experiment agree with the scale factor of the 
accelerometer for higher acceleration inputs. Application of the mass attraction principle as a 
calibration method of accelerometers for very low accelerations (for instance, in orbit) is 
proposed. 


17. KEY WORDS 


Low g accelerometer 
Mass attraction testing 
MESA accelerometer testing 


18. DISTRIBUTION STATEMENT 

Announce in STAR 


19 . SECURITY CLASSIF. (of thla report) 


20. SECURITY CLASSIF. (of thia page) 


21. NO. OF PAGES j 22. PRICE 


Unclassified 

MSFC - Form 3292 (May 1969) 


Unclassified 


46 


$3. 00 









ACKNOWLEDGMENTS 


This research was accomplished while the author held a National Research 
Council Post-Doctoral Resident Research Associateship supported by the National 
Aeronautics and Space Administration. 

For assistance in planning and executing the experiments, thanks are given 
to Mr. B. Walls who made available the equipment and supported the experiments, 
Mr. M. Vaughan who helped during the measurements and supported the test 
setup, and Mr. R. Kissel who furnished advice in computer techniques and pro- 
vided several programs. 

Special thanks are given, to Dr. W. Haeussermann, Mr, C. Mandel, and 
Dr. G. Doane III for their interest in this work and for their advice. 



TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

MASS ATTRACTION AS CALIBRATION FORCE 2 

MASS ATTRACTION OF A RIGID BODY 3 

CALIBRATION IN THE LABORATORY 7 

Experiment Setup 7 

Results and Discussion 11 

CALIBRATION IN ORBIT 24 

CONCLUSIONS 25 

APPENDIX A: DATA HANDLING FOR THE LABORATORY 
CALIBRATION OF AN ACCELEROMETER USING THE 

MASS ATTRACTION PRINCIPLE 26 

Data Analysis 27 

Data Reduction 29 

APPENDIX B: COMPUTER PROGRAM FOR DATA ANALYSIS 

AND REDUCTION 32 

APPENDIX C: MASS ATTRACTION OF THE TILT METER 42 


REFERENCES 


44 



LIST OF FIGURES 


Figure Title Page 

1. Attraction of two mass points 4 

2. Attraction of a large body 5 

3. Error for the approximation of a cube by a 

mass point 5 

4. Optimization of the attracting mass 6 

5. Test setup for the calibration of the miniature 

electrostatic accelerometer in the laboratory 9 

6. Attracting mass assembly 12 

7. Accelerometer and tilt meter data without mass 

attraction (during work time) 13 

8. Accelerometer and tilt meter data without mass 

attraction (after work time) 14 

9. Accelerometer and tilt meter data of experiment 

number 1 16 

10. Fourier analysis of the accelerometer data before 

data reduction 17 

11. Accelerometer data after data reduction 19 

12. Fourier analysis of the accelerometer data after 

data reduction 21 

13. Autocorrelation function of the accelerometer 

data 22 

14. Smoothed power spectral density function of the 

accelerometer . 23 


iv 



LIST OF FIGURES (Concluded) 


Figure Title Page 

15. Calibration of a single-axis accelerometer in 

a satellite by mass attraction 25 

C-l. Mass attraction of the mercury in the tilt 

meter 42 


v 



LIST OF TABLES 


Table Title Page 

1. Acceleration by a Lead Sphere 4 

2. Attraction of the Proof Mass by the 

Lead Body 7 

3. Acceleration of the Accelerometer Float 

by the Lead Body for Different Vertical 

Positions 8 

4. Results of the Experiments 20 


VI 



DEFINITION OF SYMBOLS 


Symbol 

g 

Mg 

ng 

X, Y, Z 
M 
m 
F 

y 

R 

A 

A (I) 

X 

r 

P 

d 

k 

E 

NA 


Definition 

cm 

Earth attraction 981. 0 — 5 - 

s 

Micro g = i0~ 6 g 
Nano g = 10 -9 g 

Cartesian coordinates of the attracting mass 
Attracting mass 

Proof mass of the accelerometer 

Force of mass attraction 

Universal gravitational constant 

Distance between two mass points 

Acceleration by mass attraction 

Variable acceleration input 

Direction of the accelerometer input axis 

Radius of a sphere 

Density 

Distance between center of proof mass m and 
surface of attracting mass M 

Length of a lead cube 

Error of replacing the cubes by mass points 
Amplitude of the acceleration input 


vii 



DEFINITION OF SYMBOLS (Continued) 


Symbol 

I 

MESA 
AX (I) 
AY (I) 
X(I) 

Y (I) 
PSD 
A (K) 
B(K) 

K 

XSM(I) 
RX (L) 
L 

F(K) 

G(K) 

GR(K) 

GO(K) 

CR(K) 

XSD 


Definition 

Measurement number 
Miniature electrostatic accelerometer 
MESA data 
Tilt meter data 

MESA data with zero mean value 
Tilt meter data with zero mean value 
Power spectral density 
Fourier coefficient 
Fourier coefficient 
Harmonic number 

Mean value of X(I) between I-JA and I+JA 
Autocorrelation function of X(I) 

Lag number 

Frequency of the Kth harmonic 

PSD amplitude of the random and periodic data 

PSD amplitude of the random data 

PSD amplitude of the periodic data 

Fourier amplitude for the Kth harmonic 

Standard deviation of X (I) 


viii 



DEFINITION OF SYMBOLS (Concluded) 


Symbol 

Definition 

AM 

Mean value of the accelerometer response to 
the input 

SC 

Scale factor of MESA 

XM, YM 

Mean value of AX (I) and AY (I) , respectively 

N 

Number of measurements 

FN 

Nyquist frequency 

H 

Sampling time 

FO 

Frequency of the acceleration input caused by 
mass attraction 

KXY(L) 

Cross-correlation function of X(I) with Y(I) 

GS(K) 

Power spectral density function 

ML 

Maximum lag number 

P 

Multiplication factor of the tilt meter data 

XR(I) 

Data X(I) after reduction 


ix 



TECHNICAL MEMORANDUM X- 53857 


ACCELEROMETER CALIBRATION 
IN THE LOW g RANGE 
BY MEANS OF MASS ATTRACTION 


SUMMARY 


Mass attraction is used as an equivalent acceleration input to calibrate 
an accelerometer. The upper limit of the acceleration by a reasonable mass 
size is 10 -9 g in orbit and 10 -7 g in the lal oratory. The calibration has been 
carried out in the laboratory for an electrostatic suspended single-axis accel- 
erometer (MESA) with a variable mass attraction. The mass attraction input 
to the accelerometer was a sine wave with the amplitude of 23 nano g. The 
response of the accelerometer to this acceleration input by mass attraction 
was obtained by data analysis and data reduction with a computer. The results 
of the experiment agree with the scale factor of the accelerometer for higher 
acceleration inputs. Application of the mass attraction principle as a calibration 
method of accelerometers for very low accelerations (for instance , in orbit) 
is proposed. 


INTRODUCTION 

Generally, the single-axis accelerometer, used for inertial navigation, 
is tilted in the earth’s gravitational field for its calibration between ±1 g and 
±i Mg, where the acting acceleration is a function of the angle between the 
sensitive accelerometer axis and the local vertical. The input axis is almost 
horizontal for low g calibration, and the accuracy of this method is limited 
by the measurement of a very small angle. The resolution of the best available 
theodolite is 0. 1 arc sec with an accuracy of 0. 2 arc sec. The angle of 0. 2 arc 
sec limits the calibration accuracy to 1. 0 /ug. This is sufficient for almost 
every inertial navigation system. An acceleration error of i. 0 jug in an inertial 
navigation system would cause a position error of 50 m after i hour. 

However, some necessary measurements in space require a calibration 
of the accelerometer far beyond this 1— /ug limit. The Apollo Application Program 
includes the measurement of the gravity gradient anomalies of the moon from a 
satellite [ 1 ] . Gravity gradient is a function of mass distribution, and every 



mass concentration near the surface of the moon will generate a deviation of 
the measured gravity gradient from the calculated gravity gradient of a homo- 
geneous sphere. These gravity gradient anomalies of the moon are expected to 
be very large because of the big mass concentrations (mascons) , which have 
already been detected by the lunar orb iter data [2], The gravity gradient of 
the moon in the proposed 55-km orbit is about 2x 10 3 Eoetvoes units. The 
required performance of the experiment is to measure 0. 5 Eoetvoes units to 
determine the magnitude and direction of the lunar gravity gradient anomalies. 

i S' 

One Eoetvoes unit is iO -9 or 1. 0i8x 10 -12 — “ — . If an accelerometer is 

sec cm 

used as a gravity gradient sensor at a 2-m distance from the mass center of 
the satellite, the required threshold is iO -10 g. 

The accelerometer must be tested and calibrated for this low g appli- 
cation. Because a range of just i0~ e g is in the state of the art, the acceler- 
ometer must be switched in the orbit to this high sensitive range after the large 
acceleration of the launch. 


An attempt has been made to calibrate an accelerometer in earth orbit 
with a centrifuge [3], (This calibration was limited to IO -6 g. ) The apparent 
acceleration is a function of the distance from the rotation center and of the 

rad 


angular velocity. The angular velocity would be as small as 10 


-3 . 


sec 


for a 


i0" 9 -g acceleration of a mass point at a 2. 5-cm distance from the rotation 
center. It would be very difficult to get such a low constant angular velocity 
in the satellite. 


This paper describes the use of the mass attraction of a rigid body for 
the calibration of an accelerometer in the range below 10 -9 g in a satellite and 
below iO -7 g in the laboratory. The size of the mass that can be handled is 
the limiting factor. 

Calibration in the laboratory with the mass attraction method has been 
carried out and is described. A proposal is made for the calibration of an 
accelerometer in a satellite. 


MASS ATTRACTION AS CALIBRATION FORCE 

Generally, a component of the mass attraction force between the 
accelerometer and the earth is used for calibration, but, for the low g range, 
the mass attraction of a smaller rigid body has some advantage because the 
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acting acceleration is a function of the size of the mass. It is very easy to 
generate a small acceleration by mass attraction, but there are two problems: 
The accelerometer must be able to sense this small acceleration, and the 
system must often separate this small input acceleration from much larger 
background disturbances. 

The lower limit of acceleration sensing is given by the threshold of the 
accelerometer and the requirements for the experiment. The upper limit is 
given by the size of the mass. Because the mass must change its position to 
generate a variable additional acceleration input, the mass should be of a 
size that can be handled. The mass used in the laboratory experiments was 
about 1000 kg; in a satellite, it might be about 1 kg. Therefore, the corre- 
sponding limits of the mass attraction are about 0. 05 Mg in the laboratory and 
0. 00 i Mg in the satellite. 


The disturbances in the laboratory are the unknown tilt angle of the 
accelerometer in the earth field, tilting of the foundation, bending of the fixture, 
and changing of the temperature; the disturbances in a satellite are mass 
attraction of some other parts of the satellite, changing of the temperature, and 
acceleration by gravity gradient of the orbited planet. Because generally 
these disturbances cannot be eliminated, they are separated by filtering as 
discussed in this paper. The input acceleration is changed with a certain 
frequency, which should not correspond with the frequency of any other 
acceleration. The input frequency of the mass attraction can be very low. The 
variation of the mass attraction is either a change in the distance between the 
accelerometer and the attracting mass or a change in the angle between the 
sensitive accelerometer axis and the vector of the mass attraction. 

The response of the accelerometer to the mass attraction force can be 
obtained by a conventional data analysis utilizing a computer. 


MASS ATTRACTION OF A RIGID BODY 

The accelerometer to be calibrated is a single-axis device; therefore, 
only the acceleration component along the input axis is of interest. The proof 
mass of the accelerometer is m and the calibrating large mass is M. The 
Cartesian coordinate system X, Y, Z has its origin at the center of m and its 
X-axis along the sensitive axis of the accelerometer (Fig. 1). 

If m and M are two mass points with the positions (0,0,0) and 
(X, Y, Z) , the attraction force F in the X direction would be [4] 
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mxMxX 


(1) 




7=6. 67xl0 -8 cm 3 gm -1 s -2 
R = (X^Y^Z 2 )^ 
m = proof mass 
M = attracting mass; 


Figure i. Attraction of two 

mass points. y is the universal gravitational con- 

stant in CGS units, and R is the distance 
between the mass points. The acceleration A of the proof mass along the 
input (X) axis is 


A = 7 M 


X 

R 3 


( 2 ) 


This equation is also correct if the calibrating mass M is an extended 
homogeneous sphere. Then R is a sum of the radius r of the sphere and 
the distance d between the center of the proof mass and the surface of the 
sphere; r changes with the mass M. The acceleration of a lead sphere on 

the X-axis with a density of p = 11. 34-^ — r in a distance of d = 10. 2 cm from 

cm 0 

the proof mass is shown in Table 1. For an acceleration of 10 -7 g the mass 
would weigh 4700 kg. This is already too much to handle easily in a laboratory. 


TABLE 1. ACCELERATION BY A LEAD SPHERE 


Mass (kg) 

Acceleration (g earth) 

Radius of the Sphere (cm) 

0. 15 

10-1° 

1.5 

4 

10" 9 

4.4 ' 

68 

10“ 8 

11. 25 

4700 

10“ 7 

46.3 

1. 4xl0 6 

i0 -6 

307 

1. 4xi0 24 

1 

3. 07x1 0 8 
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E (dimensionless) 


The mass attraction of a' rigid body with an irregular shape is given by 
a triple integral and has no general analytical solution; but if the distance 
between the proof mass and the rigid body is large in relation to the size of 
the body, then the attraction is the same as that of a sphere or a mass point with 
the mass M [5]. Therefore, the mass M was considered to be cut into 
small cubes of the dimension k with an imaginary point mass M. in the center 

of each cube (Fig. 2). The mass attraction is computed for every cube and the 
sum of all cubes is the attraction of the whole body. 



Figure 2. Attraction of a large body. 


0.1 
0.01 
0.001 
0.0001 
0.00001 

0 1 2 3 4 5 6 
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Figure 3. Error for the approxi- 
mation of a cube by a mass point. 


The quality of this approximation 
is shown in Figure 3. The error E is 
a function of the ratio of the distance 
R to the length k of the cube, 


F _ , _ A (cube) 

A (point) 

(dimensionless) . (3) 

For a 10~ 3 accuracy of the approximation 
of a cube by a mass point, the cube 
length k should be less than one- 
half the smallest distance. In the 
experiment, the distance between the 
center of the nearest cube and the 
center of the proof mass was 10. 2 cm. 

The length of the cube was 5. 1 cm. 

In the computation of the mass 
attraction, the cylindrical proof mass 
of the MESA was also considered as a 
mass point. The error of this approx- 
imation is less than 10 -3 for all cubes. 

The attracting mass was limited 
to less than 1500 kg by the available 
suspension capability. Therefore, an 
optimization was made by changing the 
shape of the mass. The mass attraction 
for every possible cube in a space of 
20k by 20k by 1 0k was calculated. The 
first 900 most attracting cube locations 
were chosen for the attracting mass of 
the experiment. Figure 4 shows the 
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Acceleration (fig) 



Figure 4. Optimization of the attracting mass. 
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relation between a sphere and an optimized body during assembly. The dis- 
tance d between the proof mass center and the surface of the attracting body 
was 10. 2 cm in this graph. 

The input acceleration is varied by changing the distance d between 
the proof mass and the surface of the attracting body (Table 2). The body 
that was used is a compound of 100 lead cubes with 10. 2-cm lengths each and 
100 lead cubes with 5. 1-cm lengths each. To get a change in the attraction 
of 70 ng, the mass must be moved about 76. 2-cm along the input axis of the 
accelerometer. 

TABLE 2. ATTRACTION OF THE PROOF MASS BY THE LEAD BODY 


Distance d (cm) 

Mass Attraction (ng) 

5. 1 

81. 8 

10. 2 

67. 2 

12.7 

60. 8 

25.4 

38. 1 

38. 1 

25. 3 

50. 8 

17. 7 

63. 5 

13. 0 

76. 2 

10. 0 


CALIBRATION IN THE LABORATORY 
Experiment Setup 

The maximum acceleration produced by a mass less than 1500 kg is 
not more than 75 ng for a nearest distance of 7. 6 cm between the center of the 
accelerometer proof mass and the surface of the attracting mass. The same 
amount of acceleration would appear by an input axis tilt of 0. 015 arc sec into 
the gravity vector. The uncontrolled movement in the laboratory of people 
stepping on the test pad causes disturbances much larger than 0. 1 Mg; there- 
fore , almost every experiment was made after working time when the test 
room was locked. But even then, the tilting of the test pad was measured 
with an accurate tilt meter during the experiments. 
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The maximum change of 75 ng could not be obtained in the laboratory 
because the necessary movement of the lead mass was impossible. The 
movement of the lead mass on the floor would cause some additional tilting of 
the test pad. Therefore , a suspension of the mass on an I-beam of the ceiling 
was considered best. But in this case, the colinear movement of the mass 
along the horizontally aligned accelerometer input axis would require too large 
a force pulling horizontally. The chosen change of the mass attraction along 
the accelerometer input axis was a movement along the vertical. 

The vertical movement changed the acceleration component along the 
input axis because the angle between the input axis and the mass attraction 
vector was varied. The lead mass was lowered and raised with an electrical 
hoist. The acceleration of the accelerometer proof mass caused by the mass 
attraction of the lead body is shown in Table 3 for different vertical positions. 
The distance between the center of the accelerometer float and the surface of 
the attracting lead body, in the zero position, is 7. 6 cm. Figure 5 is a drawing 
of the experimental setup for the calibration of the accelerometer using this 
mass attraction scheme. 

TABLE 3. ACCELERATION OF THE ACCELEROMETER FLOAT BY THE 
LEAD BODY FOR DIFFERENT VERTICAL POSITIONS 


Vertical Position 
(cm) 

Acceleration 

(ng) 

Vertical Position 
(cm) 

Acceleration 

(ng) 

0 

74. 017 

27. 9 

45. 979 

2. 5 

73. 469 

30. 5 

41. 127 

5. 1 

72. 559 

33. 1 

36. 372 

7.6 

71. 280 

35. 6 

31. 954 

10. 2 

69. 623 

38. 1 

28. 003 

12.7 

67. 571 

40.7 

24. 548 

15. 2 

65. 103 

43. 2 

21. 562 

17. 8 

62. 196 

45.7 

18. 904 

20. 3 

58. 823 

48. 2 

16.788 

22.9 

54. 970 

50.7 

14. 890 

25.4 

50. 658 

53. 3 

13. 254 


The experimental setup allowed a total vertical travel of just 30. 5 cm. 
To get the maximum change in the mass attraction, an optimization study was 
made. It showed that if the mass is moved vertically between 12.7 cm and 
43. 2 cm below the zero position, the change in acceleration is maximized. 


8 








Figure 5. Test setup for the calibration of the miniature electrostatic 
accelerometer in the laboratory. 


For the harmonic determination, the input acceleration for as few as four 
points of a sine wave of vertical position were recorded: middle, maximum, 
middle, minimum. The maximum position was at a vertical distance of 
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12. 7 cm below the zero position, the minimum position was 43. 2 cm below 
and the middle position was 28. 6 cm below. The middle position was chosen 
so that the input acceleration at that position is the average of that occurring 
at maximum and minimum positions. The succeeding run through these 
positions produces a sine-wave input A (I) with an amplitude NA 

NA * (23. 0 ± 0. 2) 10” 9 g A (I) = NA sin ( “- I) 


With a movement larger than 100 cm, the input amplitude could be as high as 
0. 037 jug with the same mass and the same components. 

To demonstrate low g accelerometer calibration, the MESA, manufac- 
tured by Bell Aerosystems, was used [6], The MESA is a single-degree-of- 
freedom accelerometer with an electrostatically suspended proof mass. This 
proof mass is a thin-walled beryllium cylinder with a flange for pickoff and 
restraint. The cylinder length is 2. 915 cm and the inner diameter is 1. 268 cm. 
The distance between the center of the float and the mounting surface of the 
accelerometer is 3. 134 to 3. 190 cm. 

The proof mass suspension force is adjustable to correspond to different 
accelerations of the environment. A pulsed force rebalance technique is 
used to constrain the proof mass along its sensitive axis. The pulse rate is 
proportional to the acceleration along the cylinder axis (input axis). The 
MESA that was used had two ranges: low g range with 10. 21 pps/Mg and 
high g range with 1. 024 pps/Mg. For the whole experiment, the MESA was 
used in the low g range. The technical data in the acceptance test of the 
MESA in the low g range are 


null stability = 0. 163 x 10“ 6 g (4-hour period) 
scale factor stability = ±0. 01758% (± 0. 1% - 4 hours) 

The mounting of the MESA is a clamp type around the case to get a 
small distance between the MESA and the lead surface. 


A tilt meter made by Ideal Aerosmith was set on the same test pad to 
measure tilts of the foundation [7]. The surfaces of two interconnected 
mercury pools 1 m apart serve together with two rigidly mounted plates as 
two capacitors. Any tilting will change the capacities , which can be expressed 
in tilt angles. The tilt meter was used in the high range. The output of the 
null meter was calibrated so that a tilt of 0. 0067 arc sec produced a voltage 


10 



of 0. IV. This voltage was fed to an integrating digital voltmeter. Both the 
accelerometer data and the tilt meter data were automatically punched on a 
paper tape for computer use. 

The attracting mass, 100 lead cubes of 10. 2 cm lengths and 100 lead 
cubes of 5. i-cm lengths, was assembled in the optimized shape (Fig. 6). The 
free spaces in the basket were filled with the same size of wooden cubes. 

The time for every measurement was H=100 seconds. This large 
measuring time was necessary to get a minimum number of impulses required 
for statistical data handling. During the measurement, the mass remained in 
position. At the end of one measurement, the mass was raised or lowered 
to the next position. The time between the measurements was always 20 
seconds. In every series, at least 100 measurements corresponding to 25 
cycles of a sine wave were made; 


A (I) = (23. 0 ± 0. 2) sin ( “ I) ng 
I = measurement number 


The input acceleration is given in ng = 0. 981 x 10 6 -—jf- . The experiment 
had to show the response of the accelerometer to this small input acceleration. 

The constant acceleration (or bias) input to the accelerometer, perhaps 
caused by tilt, was about 17 Mg. The setup is very sensitive to temperature 
changes; therefore, the temperature was controlled in case II (Fig. 5) and 
was constant to less than 0. 01° C for experiment number 1 (less than 0. 1° C 
for other runs). 


Results and Discussion 

The experiment had to show how to measure the small acceleration caused 
by mass attraction and how to define the scale factor of the accelerometer for 
low acceleration. Because the acceleration input was a relative change, the 
constant part of the data was removed from the accelerometer data AX (I) and 
the tilt meter data AY (I). The accelerometer data X(I) and the tilt meter 
Y (I) have zero mean value. These data are used for the computer analysis. 


The random noise is much higher during work time than after work 
time. Figure 7 shows the random noise during work time. There was no 
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Figure 6. Attracting mass assembly. 

movement of the lead mass in this experiment (run number 6). The data were 
collected during work time every 2 minutes with a 100-second measurement 
time as in the later runs with mass movement. Run number 5 is the same as 
run number 6 just after work time (Fig. 8). The standard deviation from an 
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Acceleration (r*g) 





adjusted curve for the run during work time is twice as much as for the run 
taken after work time (all the experiments were made after work time except 
run number 4). 

The response of the accelerometer and the tilt meter to the mass 
attraction is overlapped by stochastic and other periodic acceleration. But 
already in the raw data, the movement of the mass is clear (Fig. 9). The 

data X (I) and Y (I) are shown in Figure 9 for run number i after subtracting 
the mean value. The time scale includes the necessary 20 seconds for moving 
the mass. The tilt meter data were multiplied by a scale factor to get radians 
that correspond to acceleration in g. The dimensions for the accelerometer 
data were originally impulses per second. But to give a better impression 
in all the graphs, the scale factor of the accelerometer data sheet is used to 
express the data in pg. For every experiment the first position of the lead 
mass was the middle position followed by the upper position. 


The tilt meter data also show a very small modulation with the frequency 
of the mass movement. However, it is much less than the corresponding 
response of the accelerometer. If the accelerometer data were the result of 
some tilting corresponding to the mass movement, the tilt meter data would 
show this. The small modulation of the tilt meter data is caused by the mass 
attraction and is computed in Appendix C. The tilt meter is considered as a 
supporting device and the analysis is shown only for the accelerometer data. 

The accelerometer data are digital, and the analysis is made with a digital 
computer. A statistical analysis (autocorrelation function and power spectral 
density) is used together with a Fourier analysis to compute the response of 
the accelerometer to the small mass attraction input. The detailed analysis 
and data reduction techniques are given in Appendix A , and a computer program 
is given in Appendix B. 

The acceleration input is a sine wave with a frequency of 0. 0025 Hz. 

This input frequency corresponds for 100 measurements to a harmonic number 
of K=25. The unforced response evaluated for a sine wave with a frequency 
of 0. 0025 Hz is the Fourier coefficient B(25). For a sample of 100 data X(I) 

100 27T 

B(25) = 0. 02 2 X(I) sin (— x I) 

1=1 4 


where X(I) is the accelerometer data with zero mean value and I is the sample 
number. The accelerometer data X(I) in Figure 9 show the shape of a saw 
tooth with the frequency of 0.0025 Hz. Therefore, the Fourier analysis of the 
data without a data reduction has large amplitudes with low frequencies 
(Fig. 10). The Fourier coefficient B (25) is not just the response to the input 
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Figure 10. Fourier analysis of the accelerometer data before data reduction. 


sine wave because a certain part of the 25th harmonic is produced by the saw- 
toothed shape of the overall data. Therefore, several types of data reduction 
were tried (Appendix A ) ; two of which were used for the calculation of the 
following scale factor: 

1. Subtraction of XSM(I) from every X(I). XSM(I) is the mean value 
of the data X(I) between I-JA and I+JA. For JA=3, the summation was made 
over 7 points plus two edge points with a weight factor of 0. 5. 

2. For the error calculation, the standard deviation was computed after 
subtraction of the 20 highest and 20 lowest harmonics. The accelerometer 
data X(I) are shown in Figure ii after the subtraction of the adjusted curve 
XSM(I) from X(I) with JA=3. The amplitude of the input acceleration was 

23 ng at a frequency of 0. 0025 Hz. 

The Fourier analysis CK(K) gives the amplitude as a function of the 
frequency. In Figure 12, the amplitudes of two Fourier analyses are. shown. 
The one amplitude distribution is from a run without a mass movement. This 
experiment is run number 5 in Table 4. The amplitude at 0. 0025 Hz is almost 
the same as the mean value of the amplitudes at the frequencies from 0. 0020 
to 0. 0030 Hz. The other curve in Figure 12 is the amplitude distribution of 
a run with a mass attraction input at a frequency of 0. 0025 Hz. This run is 
number 1 in Table 4. The Fourier coefficient B(K) is 27. 80x 10 -2 pps with 

a frequency of 0. 0025 Hz. The mean value of the amplitudes of the harmonics 
from K=2i to 24 and K=26 to 29 is as low as 2. 17 x Hr^ps. The amplitude at 
K=25 in run number 1 is supposed to be a summation of the response to the 
sine-wave input and of some additional random noise. Because the bandwidth 
is limited, the subtraction of the mean random noise level from the amplitude 
at K=25 is justified. Therefore, the response of the accelerometer to the sine- 
wave input at 0. 0025-Hz frequency is 25. 63 x lO-^ps. This value is very 
close to the computed one with the scale factor for higher acceleration, 

23x 10 -2 pps. 

The statistical analysis shows almost the same results as the Fourier 
analysis. The autocorrelation function RX(L) of the accelerometer data 
X(I) shows the periodic response to the input acceleration (Fig. 13). From 
the autocorrelation function, the power spectral density is computed. The 
amplitude G(K) of a periodic function is given by 


G(K) = \ r GX(K) XB , 

where GX(K) is the power spectral density function and B is the bandwidth. 
The distribution of the amplitude G(K) of a smoothed power spectral density 
is shown as a function of the frequency in Figure 14 for run number 1. The 
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Figure 11. Accelerometer data after data reduction. 


TABLE 4. RESULTS OF THE EXPERIMENTS 


Experi- 

ment 

No. 

NA 

Input 

Accel 

(ng) 

Accelerometer Response 

Random Noise Part 

Signal Response 

Standard 
Deviation 
10 -2 pps 

G(25) 
10 -2 pps 

B (25) 
10 -2 pps 

GR (25) 
10" 2 pps 

CR(25) 
10 -2 pps 

GO 

10 -2 pps 

BO 

10 -2 pps 

1 

23 

28. 07 

27. 80 

4. 70 

2. 17 

23. 27 

25. 63 

5. 59 

2 

23 

26. 56 

26. 80 

3. 87 

3. 41 

22. 69 

23. 39 

6. 77 

3 

23 

29. 72 

29. 64 

5. 02 

5. 57 

24.70 

24. 07 

11. 90 


Mean value for No. 1, 2, 3 


24 ± 2 


23 


24. 78 


25. 62 


5. 85 


5. 30 


18. 93 


20. 32 


14. 06 


5 

6 


0 

0 


2 . 02 

3. 14 


1. 12 

2. 08 


2. 14 

3. 94 


2 . 01 

4. 13 


-1. 2 
-0. 80 


-0. 89 
-2. 05 


4. 90 

10. 12 






RXX 





amplitude at the frequency 0. 0025 Hz is G(25)=28. 07 x i0 _2 pps. The mean 
value of the random noise for the frequency FO is GR(25)=4. 70xi0~ 2 pps and 
the response to the sine wave at the frequency of 0. 0025 Hz is therefore 
23. 27 x 10“ 2 pps. 

For the error calculation, the 20 higher and 20 lower harmonics of the 
Fourier analysis are subtracted. This does not change the amplitude of the 
25th harmonic. After this filtering, the standard deviation from the ideal 
response sine wave with the amplitude B(25) is calculated. This is 
XSD = 5. 59 x 10~ 2 pps for run number 1, which was made after work time with 
a very low noise level. Runs 2 and 3 were made under almost the same con- 
ditions, except that the noise level was higher because of the air-conditioning 
equipment. Run number 4 was made during working hours; because of the 
high noise level, the data could not be used for scale factor calculations. 

The mean value of the response to the mass attraction input is for the 
first three numbers, 


AM = (24 ± 2) x 10 ^ps . 


This gives a scale factor SC=(i. 04±0. 09)xi0 7 pps/g for the input amplitude 
of 23. 0 ng. This result is very close to the scale factor of the data sheet in 
the higher acceleration with 1. 02xi0 ? pps/g. 


The size of the error for the scale factor could be reduced by an improved 
setup incorporating better noise isolation and automatic and larger movement 
of the attracting lead mass. For larger movement of the mass, the amplitude 
of the input sine wave could be twice as high, and the calibration would be more 
accurate. The experimental setup was more or less improvised. The aim of 
the experiments was to prove that the calibration with a mass attraction in the 
low g range of 10 -8 g is possible with good confidence and repeatability. The 
application of these experiments to the calibration of an accelerometer in a 
satellite for the measurement of the gravity gradient anomalies of the moon is 
proposed. 


CALIBRATION IN ORBIT 


The accelerometer gets a large acceleration during the boosting phase 
of the rocket. Therefore, the state of the art of the accelerometer requires 
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a switching to the highly sensitive mode after the start of free flight. After 
this point, there are no large accelerations and the mass center of the sat- 
ellite is almost under 0-g conditions. The threshold of the accelerometer can 
be set very low , but there is no way to calibrate the accelerometer in this high 
sensitive range with the generally used terrestrial means. There are many 
different unknown acceleration inputs to the accelerometer in a satellite; that 
is, gravity gradient acceleration of the orbited planet, mass attraction by the 
satellite, low thrust impulses, etc. This paper proposes the calibration of 
the accelerometer in the satellite by using a well known but variable mass 
attraction scheme. A reasonable size of the mass for the calibration in a 
satellite might be 4 kg, which sets the upper limit of the acceleration to 10~ 9 g 
(Table 1). The required threshold of an accelerometer in the lunar gravity 
gradient anomalies measurement is lower than iO~ 10 g. The acceleration 
input caused by mass attraction can be varied by changing the distance or the 
component of mass attraction along the accelerometer input axis. The changing 
of the component has some advantages. The accelerometer is fixed in relation 
to the satellite. The mass can be turned around the accelerometer continuously 
with a small angular velocity. The input acceleration is a cosine function for 
the arrangement in Figure 15. The argument of the cosine is the angle between 
the input axis and the line from the center of the proof mass to the center of 
the calibrating mass. For easy data handling, the rotation axis of the cali- 
brating mass should go through the center of the accelerometer proof mass. 

If there is no other time variable input, the response of the accelerometer is 
a cosine function with an additional constant caused by the constant mass 
attraction input of the satellite. With a small computation, the scale factor, 
nonlinearity, and constant bias of the accelerometer can be calculated. The 



Figure 15. Calibration of a single- 
axis accelerometer in a satellite by 
mass attraction. 

calibration scheme for accelerations 1 


movement of the mass around the 
accelerometer can be a continuous 
one or in steps to facilitate digital 
analysis. The calibration does not 
interfere with the experiment and can 
continue through the whole measure- 
ment. 

CONCLUSIONS 

The calibration of an acceler- 
ometer in the laboratory by the mass 
attraction method may be of interest 
in the future , when the low acceler- 
ometer thresholds are required. In 
orbit, the mass attraction scheme is 
the only presently known practical 
>w 10" 9 g. 
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APPENDIX A 


DATA HANDLING FOR THE LABORATORY 
CALIBRATION OF AN ACCELEROMETER USING THE 
MASS ATTRACTION PRINCIPLE 


The measurements of every experiment are numbered from 1=1 to N. 

The input acceleration A was a sine wave with KO»4 positions per wave length 
(cycle) 


A (I) = NA sin ““ (I-i) 

I = 1,2 N 

NA = 23. 0 x iO -9 g 

N = 100 = number of measurements 

KO = 4 = number of measurements per cycle. 

There are two equal sets of N data; the MESA data, AX (I) (an integer number 
of impulses per H = 100-second measurement time) and the tilt meter data AY(I) 
(an integer number from the integrating digital volt meter). The mean values 
of the original data are 

1 N 

XM = — S AX(I) 

W 1=1 

1 N 

YM = — 2 AY (I) 

N 1=1 

for I = 1, 2, . . . , N 
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The data are transformed to data with zero mean value 


X(I) = AX (I) - XM 
Y (I) = AY (I) - YM 

These data, X(I) and Y(I), are used for the analysis of the MESA data, AX(I). 


Data Analysis 


The highest frequency in the data should be lower than the Nyquist 


frequency (FN) , which is half the sampling rate. 


Thus 


FN = 


2H 


and 


H = sampling time. The sampling time H was for every measurement 100 
seconds, giving a Nyquist frequency of FN = 0. 005 Hz. The input frequency 
of the mass movement is FO = 0. 0025 Hz. 


The statistical analysis uses the autocorrelation function and the power 
spectral density function. The autocorrelation function of X(I) is given by 


BX(L) = 


1 

N-L 


N-L 

2 X(I) x X(I+L) 
I = i , 


5 


when L = lag number 1 to N/2 and N = number of samples = 100. The normal- 
ized autocorrelation function RXO(L) is obtained by dividing RX(L) by RX(0). 

The cross-correlation function of the MESA data, X(I) , with the tilt 
meter data, Y(I), is 

RXY(L) = — — 

' N- L 

for L = 1, 2, . . 


N-L 

2 X (I ) x Y(I+L) 
I=i 

N 
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A raw estimate of the power spectral density is [8, 9] , 


GX(K) = 2Hx 


7rLK K 

RX (0)+2x 2 RX(L)xcos ( ■■■ )+(-l) xRX(ML) , 

L+l Mij 


where 

K = harmonic number 
H = sampling interval 
ML = maximum lag number 
L = lag number. 

The frequency for the PSD is 
TC 

=— 

v ' 2xMLxH 

and the bandwidth is 

B _ — i 

MLxH 

2 

GX (K) is given in units /B. 

The amplitude versus frequency is, for a periodic function, 

G(K) * N GX (K)xB 

For the analysis, the nonsmoothed PSD is used. 

The Fourier analysis is also used to calculate the amplitude of peri- 
odic functions. The Fourier coefficients for the MESA data, X(I), are 
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K = i 2 il 

ix 1 j " » • » • > 2 

K= 1,2 ~ 

A (0) = 0 

where 

N = number of data (even) 

K = harmonic number. 

The frequency of the Fourier components is 

F (K) = — — = Kxl0“ 4 Hz 
NxH 


A(K) 


_ 2 _ 

N 


N 

2 X(I) 
1=1 


cos 


2ttK 

N 


B(K) 


_ 2 _ 

N 


N 

2 X(I) sin 
1=1 


27rK 

N 


Because the maximum lag number in the statistical analysis was chosen to be 
N 

ML =— , the frequency of the Fourier analysis agrees with the frequency 
of the PSD for the same K. 


Data Reduction 

The amplitudes of the lower harmonics in the data analysis are larger 
than the observed amplitude at the frequency FO of the acceleration input. To 
unmask the amplitude at frequency FO, four types of data reduction are 
available in the computer program (Appendix B). First, the trace of the tilt 
meter is subtracted from the MESA data after the adjustment by a factor P. 
The factor P is calculated with the least mean square fit method [10] 


N 

S X(I)-PxY(I) 
1=1 L 


= Flin 
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SX(I)xY(I) 
2 Y(I)xY(I) 


The reduced data are 


XR(I) = X(I)-PxY(I) 
for 1=1,2, . . . ,N. 


The second type of reduction is a filtering by subtraction of the lower and 
higher harmonics 


10 

XR(I) =X(I) - 2 

K=1 



+ B(K) sin 


2ttKI " 
N . 


50 r 
-2 
K=40 


A (K) cos 


2ttKI 

N 


+ B(K) sin 


2ttKI 

N 


A (K) and B(K) are the Fourier coefficients computed in the previous data 
analysis. 


The third reduction is the most effective one to filter out the sine wave 
with the wavelength KOxH. From every data point X(I) , a mean value of the 
S closest points is subtracted to get the reduced data XR(I). 


XR(I) = X(I)- 


2S 


I+S-i 

0. 5XX(I-S)+ 2 

J=I-S+1 


X(J)+0. 5XX (I+S) 



2S must be a multiple of the sample number KO per wavelength to avoid 
changing the amplitude at the frequency FO. 

In several experiments, the data X(I) build almost a straight line 
that is modulated by a sine wave with the frequency FO. Therefore, the sub- 
traction of the adjusted straight line would reduce the data, but it is not 
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sufficient in every case. After reduction the data again go through the analysis. 
The standard deviation XSD of the reduced and filtered data XR(I) from the 
idle or unforced response to the sine wave contains the random noise and is 
considered as a good approximation of the error, 

XSD =| ~ s[xR(I) - B (25) sin ( ^ I)] 2 J 

RX(I) = reduced and filtered data 

B(25) = Fourier coefficient of the 25th harmonic. 

The XSD with the amplitude G(25) of the power spectral density gives the 
same value as with the Fourier coefficient B(25). 
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APPENDIX B 


COMPUTER PROGRAM FOR 
DATA ANALYSIS AND REDUCTION 


This program is written in FORTRAN II for the GE235 digital computer. 
The arrangement of the program and the data follows. 

1. SLEM (computer control cards) 

2. Main program 

3. Subroutine plot 

4. Data card i/2 

5. For data input 0 Format (14) 

For program test 1 Format (14) 

6. Experiment card N, NH, ML, NA, MN, MAXL, 


MAXK, KO, MON, MD 


Format (10(14) ) 


7. MESA data AX (I) 


Format (5(3X16, IX)) 


8. Tilt meter data AY (I) 


Format (5(2X18, IX)) 


9. Sample number with KM (I) 

wrong data 


Format (14) 


10. Data reduction IC 


MA or JA 


Format (14) 
Format (14) 


11. End 


-1 


Format (14) 


The subroutine plot is for plotting the data on the initial printout. 


PROGRAM LISTING FOR DATA ANALYSIS & REDUCTION 


C CALIBRATION OF THE MESA 

XOMMHN Y(100 J »XSM tlOO I .RYYl&O I »Q(50 L 

l^rerr7RrT3TT7imT51t7RYX [ 311 , GY [gn~r — 

2RXY0I50 ) . AX[600 J.KMI201 * AYIoOO] ,GXOI50],FC50 J, 

3 K L ( 6 n 0 ] iGXXt50)*AK[60],BKt‘>C! - - 

DIMENSION X[10nj,RXXIlO0],OxU00],CK(100]»GXAll0O J 

PI S3”, 141592634 

P12 = 6, 283185308 

TT'ROD ~YVT','<0 

201 FORMAT [ I 4 I 

if Tim ?i7,?rr,2-(r7 


C TEST OF THF PROGRAM 


202 READ 203»N.NH,MAXL»MAXK,LXS»KXS*LXC»KXC,LY,KY,KO 

203 FORMAT IT 14) ,3M8, Him] - 

KM f 1 ) sN + 2 

NA = 100 

ZLXS=LXS 

TKXSsTXS 

ZLXC=LXC 

— ZKXC _ KXj , : - 

Z L Y = L Y 

-?ky=ky - 

PRINT 6,LXS,KXS,LXC.KXC 

5 FURM ATT? 1TT1 TEST XTTT - I BTT5TT — SrRT7*T>T* ZIV 1 4 

16rt J + I 8 * 15H COS 1 2*P I *ZL/ 14, 3H I ] 

PRINT 18 »LY,KY 

18 FORMAT 1 22H YU) * 18, 

TT5ff“ _ CO'ST2*T J T*71.7 T4T3H ] 1 

DO 204 L = l> ML 

ZUilT-T ; 

KL(L)=ZLXS*S!NF(Pi*ZL/ZKXS)+ZLXC*COSF[Pl*ZL/ZKXC) 

~2 axtltvklt lt 

DO 205 L = l» ML 

ZUSI-1 - — 

KLtL)=ZLY*C0SFtPI*7L/ZKYJ 

~2D5~" AT rCT=KTTLT 

GO TO 206 


V READING AND "ARRANGING OF THE DATA 

210 READ I7"NTNH,ML,NA,mN,MAXL,MaXK,K0,M0N,MD 
1 FORMAT [10114] ] 

“ READ 102, [ K L [ L T TL'al.Mt - ! 

102 FORMAT [ 5 [ 3X16,1X11 

DO 71 L=1,ML 

71 AX(L)sKLtL) 

READ 103, [KLTkJ,l»l,MLJ 

103 FORMAT [ 5 [ 2X18,1X11 

DO 73 L = 1 » Ml 

73 A Y t L 1 *KL t L 1 


PROGRAM LISTING (Continued) 


C CALIBRATION 'OF f THgi M6SA 


READ 201 > (KM 1 1 } t til* MN’} 
206 ZNsN 
HsNH 

ZKO*KO 

ZML=ML 

zwmrA ^ 

ZMN*MN 
I Col 

PRINT 26 

-26 roRFsi mm 

PRINT 106 # .[AX[Ll#Lal»MLJ 

106 FORMAT [5[F12,1J I 

PRINT 106, t A Y t L J *L b 1»ML1 
PRINT^Ol. t KH [ L 3 » Lsl * MN ] 

DO 11 I si , MN 

K = KH ( I ) 

MsmL- I 

DO 11 LsK.M 

AXTTT=AX tL + l ) 

11 A Y [ L ] = A Y ( L+l I 


MeMAXL 

zitsn 

ZMAXKsMaXK 

rHAnrtsHxxr, 

A = 0 

BT ='0 — 

DO 17 1=1. N 

xt n = Axm 

A = A + AX I I I 

- IT -B=e+TYTT I ‘ : 

X H = A / Z N 

YH=B77N " 

D 04 1 = 1 , N 

A X t I T = A'XT I I -XM ~ 

4 Y (I ] = A Y [ II -YM 

216 READ 201, IC 

' GO TC t 16»119,T29Vl3T,T49,13r, m 


C DATA ANALYSIS 


16 lD=in+l 

18=10 O^ID+IC 

■ ASTI 
8 = 0 

DO 2 0 TsI'.N 

xm = x m-x‘i 
a=a+x m*xt n 
20 8 = R + Y(IJ*Ym 
RXO=AV7N 
R Y O = p / 7 N 



PROGRAM LISTING (Continued) 


c calibration or the mesa 

XSI) = S3PTF t ABSF C R X 0 J 1 

ysd=s3ktf c abspiryoi i ~ — 

PRINT F-.JH.MON.MD 

R FORMAT t 4 H 1 CAL I BPAT ION OF' THE MESA ' 14, 

15X14, XH/J2//I 

PRINT 7 4 » t X t I J # Y t IT » 1 s 1 » N ) 

74 FORMAT t2lE12,4]3 

CALL Pi 0 T [ X , N 1 - 


c autocorrelation OF XIII 
00 24 l = 1,M 
K a N - 1.. 

ZK = K 

DO 21 1=1, K 

I L = I + L ' 

21 A = A + X U ] *Xl III 

— 27 RXIL 1 = A/ZK 

24 RXX(L3=RX(L]/RX0 


C AUTOCORRELATION OF Yt I I 

pg -j 4 l-i 

KaN-L 

7KgK 

A = 0 

DO 31 1=1, K - 

I L= I +L 

T1 _ A s X+Y rrprmT — 

32 RYILI=A/ZK 

34 - RY Y TU = RYCL J7RY0 

B=RX0** , 5 

TTaTTYn**;? 

D = B * C 


C CRTTSE'C 0 R R E rATyOTr OF^XTH AWTTn T 

W'RprTTaTT'M' 7 

Qt'LJaL 

(TsTFi: : 

As 0 

l ^- K 

DO 41 1=1, K 

1 L=T+L" 

41 A8A + xm*YCIU 
RTTYTOW7K~ 

44 RXY0IL)=RXY[L3/D 

; ’ PRINT 8, IB. MON, MU ' 

PRINT 81 

81 FCRRTT [TTQ H L EH RY RTY ROOT 

1 RXY RXYO //I 



PROGRAM LISTING (Continued) 


C CALIBRATION OF THE MESA 


RC = 0 

U'J 33 lsl.W ; 

33 RC=RC+X( 11‘VIII 

rcotrctd ' 

A = 0 

8 -=X"."n 

PRINT 30. A,RX0»RY0,RC,R,B,RC0 

-m— formth — r^OTTrsTErzrr J » 2 x s t f i n . 5T1 

DO 84 l. = l.M 

84 prin t 3n. — arc i ,~ r x r t tv rt t t ) , r xrtr i . r xy tun r y y i l i » r xrortrr 

PRINT 112 


— — — -TtT 

FORMAT t XI 
PRINT 2.XH.YM 



~T 

“FORMSTteM— XM = El 2 3 > 
PRINT 12.XSD.YSD 

6 H 

TK - B1 2 VA J 

tT 

FORMAT IbH XSD= ElTTTr 
CALL PLOTIRXX.M] 

6 H 

YSD= E12.4///T 


r p ot e r s p e c t rum ar ott n 


PRINT K.IB.MON.MD ' 

PRINT 91 

91 F ORM A T 1 6 4 H K FREQUENCY RAW PSD AMPLTT SMOOTH PSD 

1 AHp tl T // 3 

MK = M-l 

Tff= 1 »r — 

BF=1 , 0/ I ZM*H I 

DO“5nC=TrtT 

Q I K I =K 

—— 

F[K]=0.5*ZK/tZM*Hl 

— Ar^O ’ ‘ 

B = 0 

DO~ 53 L = r. HK 

ZL S L 

~ 53 AsA + RXH)* C OS r [ Pl»ZL»ZK/ 7fn ' 

GXtK]g2.0«H«IRXO+ 2.0*A+RXtM)»l-1.0]»*Kl 

AB= 0 

rx=n 

DO 109 L = 1*MK 

BX=0.5*[AB + GXllIl 

A B = 2 , ll*H* I R X0* 2 ,0*AB+RX[M)1 ' 

A XOs[ABSF(BX*BF * 2 , 011**0.5 

S'sti : 

SsSQRTF [ ABSF tAB»BFI 1 

PRINT 19.A. A. AB.S.BX. AXO ~~ 

19 FORMAT ( F6 . 1 # 5 I E12 . 4 ) 1 

GXXIli= 0,5*GXI11+0.25*IAB+GX[2I] 

DO 311 K=2,MK 

KR = K _ 1 

KP=K+1 


3 $ 



PROGRAM LISTING (Continued) 


C CALIBRATION OF THE MESA 


311 GXXlK) = ri.25*GX(KRl + r;.S*GX!K I +0 . 25 *0X f Kf I 

gxx tH)8D. - s«GXiMKi Tr;-9ytrmr}- — 

D O 312 K=i,M 
DO 93 Kb 1 , M 

93 PRINT 19 > — TJTK TTF I K 1TT3TTK T i 3 X A f KT, GXXT KTf~G XOTK T 

PRINT 112 

PRI N T 3 , B F : - 

3 FORMAT t 13H BANDWIDTH = E12.4, 5H HZ //) 

CALL PLOTIGXAyWI " ~ - 


C FOURIER ANALYSIS 


DO 403 K = 1 » 50 

ZK=K - — 

FIK]=ZK/IH*ZNJ 

— - — — — 

B = 0 

DO 4t)2 I=1»N ~ " 

Z!=I-1 

Ao A + X [ I ]*C0SFIPI2*ZI*7.n/ANl 

402 B=B+X(I ]*SINF[PI2*ZI*ZK/ZN) 

" A K [ K I = 2 , 0 *7T7TN - 

BK[KJ=2.0*B/ZN 

403 CKIKTsTAKIK J**2, 0+BKIK1** 2TTTI**0'.5 

PRINT 8, I8.MQN.MD 

PRINT 404 

404 FORMAT 1 52H K FREQUENCY AK BK CK 

DO 405 K = 1 » 50 r ; - 

405 PRINT 406,0 I K ) , F [ K ) » AK I K ) , BK I K I , CK l K ] 

406 FORMATTFTTiT, 4 [ E12 , 4 IT - 

PRINT 112 

PRTNT23 

23 FORMAT I 85H N H ML INPUT ACC 

“ 1NR OF ERRORS ‘Ma^TT TOOT! 

PRINT 105, ZN, H,ZML.ZNA,ZMN,ZMAXL,”MAXK 

105 format i7[Fi2,in 

CALL PLOT (CK* 50 I 


TT STANDARD DEVIATION OF X(l] TO "XINP 
AbO 

XINPbZNA*10,0**I'9.01 
SKIbBK 1 25 J /X I NP 
DO 127 1=1, N 

ZI , 1 . 1 

U sXm«BK(25J*S!NF[l,570796327*ZI I 

1Tr -£-4A“+U*U 

SD1=SQRTF[ A/ZN J J 
ESK1bSD 1*SK1/BKT25 J 
MSsM/2 

SK^GXATM^l/XTNP 

AbO 
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PROGRAM LISTING (Continued) 


C CALIBRATION of the mesa 


D0128 I=1.N 

z i = i - a 

U = X ( 1 ) -GxA [MSI *s INF 11.570 7 96327*ZI I 
12R A S A +1 > * U 

SD2 = SQR TF l A/Zn] 

~ ESK2=SD2*SK2/GXA [MS I 
PRINT fi , Ul, MON, MO 
FRIKiT 121, XI MR 

121 FORMAT 1 25 H INPUT ACCFl.FRAT ION A - E12-4. 

122H SINf?*PI+II-l 1/ 4 ] G >n 
PRINT 12?»BKC25J,SP1 

"122 FORMAT 1 25H OUTPUT f FOURTFR I B =[E12.4»3H + F12,4, 
1 2 8 H ] s I N ( I * P I * [ 1-13/ A 3 IMP/SEC //I 
- PRINT 12 3 V G X A [ MS) , S D 2 

123 FORMATI25H OUTPUT [PSL’I C s[fcl2.4.3H + £12,4, 

1-1]/ A] IMR/SEC //I 

PRINT 124 

-124 FHRMATT34H SCALEFACTDp FOR VERY LOW 0- INPUT //] 

PRINT 125 , SKI , ESKl 

125 FORMAT! 21 H Sc"! FOUR) = E12.4.3H + F12.4, 

114H I MP/SEC/G //I 

PRINT 126.SK2.ESK2 

126 F0RMATI21H SCTRSDI = E12.4.3M + F12.4, 

114H IMP/SEC/G //I 

GO TO 216 


— cr da nr reduct ion 

c »*** * (■ ******* * 


— C tTFA ST ' ME I'M " SQUARE FIT WITH TILT DATA 

119 , Aan • 

B = 0 

— 00 9 I=1»N 

8 = B + A Y t I I * A Y t I 1 

9 ‘ A s'A + A X T I 1 * A Y t ] ) — 

AsA/B 

g-fj 

DO If) 1=1, N 

xm=Axm-A*AYm 

10 8sB+X[ I ] 

X M = B / Z N - 

GOTO 16 


129 A = 0 

DO 400 I=l»9 

THO A* A + A X f 1 1 

XA1 = A/9 , 0 

NU = N»B 

AbO 
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PROGRAM LISTING (Continued) 


C CALIBRATION OF THE MESA 


DO 401 I=NL,N 

401 A B A+AX 1 1 1 

X A2 s A/9 , 0 

Jfjj-Q 

DO 413 I *1» N 

zm : 

X m=Axm-XAl-tXA2-XAl]*IZI-5.0JV91,0 

413 A-A + X [ I ] ; 

XM= A/ZN 

GO TO 16 


C SUBTRACTION OF THE LOWER AND HIGHER HARMONICS 


139 A = 0 

READ 201. MA 

MB=50-MA 

0TT4Tr"T=TTN ‘ 

ZIsI-1 

- FT='0 

DO 4 11 K = 1 » M A 

411 B=8+AKlKI*C0SFtPl2*ZI*ZK/2NJ+BKIK)*SINF[PI2*ZI*ZK/ZNJ 

DO 412 KSW/50' 

ZK = K 

412 B=FT+AKTK)»COSF t PT2* ZT *7K77HT +BKTKTTS1 iqTt PT2 *71 *7K7ZN J 

XC1I= XII1-B 

410" AsA + xm 

XM=A/ZN 

GOTO 16 “ " 


| SUBTRACTION OF AN ADJUSTFD CURVE 


149 READ 201, JA 

Z J A = ( J A + 1 1 * 2 — 

JB = J A + 2 

J C = 1 5 0 - JA-i ~ - 

J0=100-JA 
DO 701 I =JB , JC 
MAsl-JA 

mb=i+ja 

MAI sma-1 

"'H'BIsmB + 1 ‘ 

A = 0 

00 7 5 0 K s M A , H B 
700 A = A + AX t K ) 

70 j XSM ( I ] = [ A + 0 , 5* I AX I NAl'J + "AXTMBTJTT /Z'UA 
JA1=JA+1 

DO 703 1-3, JA1 — 

LK = 2* [ I -1 I 
ZLKs'UK 
LK1 = L < + l 

' " ' ' A = 0 

OO 702 K=2,LK 
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PROGRAM LISTING (Continued) 


~ v 

U v ' CAL I RRA r 1 0N.IOF TM| m£SA 4j 


7 02 A=A^AX IK 1 *• : • > -v 

' r •' * 

- mmk. 

70 3 X S M [ I j = [A + U # 5 *i-'AX[ 11^ AXU-fcUuXl J/Zl K 

XSM 1 1 3 = 1 5* C AXt 1 3 ^ > AXf2 H 

XSMm = l ,5*AXliJ *AKTZ J + 

DO 705 I = JD# 98 : * . 

'*?*%&. ** - 

LK=2* 1 -9V ., . ' .. ~ ;— 

ZLK»2*UOO-IJ 

■ ■ -^Sr 

. tm& • 

lki=l,:k-i ■ 

A s 0 


7 DO 704 K s PC7*r ~ • 

704 ASA + AXtKl ■ .* 


705 X 5 M l I J = lA+o*lAXl IB UT ♦ .TXTUKTTT J X/.LK 

XSM [ 99 1 3 . 5* [ AXt 99 ] f , D* t AX t 98 1 tlO 0 ) ] 1 

AnO 

AXIL r? ittt 

DO 7 0 & 1 3 1 f H 

,;■ ,, . x ( u s ax ( I j - xsmt f } ;v • ■; 

• c \. f4 ^ 

:£-m 

“ : ■ J 7 06. A=^ + )t[ I'3 

. ,, 1 , ; xfi^AxzKi ■■ ■ - f ■ 


• " . Ou rq 16 • * l 

*? A to * i* 0 , ,, . • , f 

. • .. ■ ■ 'v ■ .#■ ■■■■ j/[. 

• 217' ST8P '¥ ’ ~~ r , . - . 

ffNO- 7 ;• . ;• -• , ' - 



PROGRAM LISTING (Concluded) 


SUBROUTINE PLOT [y,N] 

DIMENSION Y ( 100 I, K 1 120)# THIN 111 I 
YMINtiJsYtlJ 
YMAX = Y 1 1 J 
00 1 I =i* N 

YMINIlJsMINlF[YMIN[l),Y[I}j 

1 YMAX = MAXlFt YMAX.Yl j] ) 

DO 2 1=2.11 

2 YMINIIJs I YMAX-Y.MINti ) J/IR.+YMINI 1-1 ) 
PRINT 3 

3 FORMAT [1H1///J 
PRINT4. (VMINUJ# Isl.llJ 

4 FORMAT ( 7 XI 1 1 1 XE9 , 2 I ) 

P'R I NT5 

5 FORMAT [14X101HJ--. -I I 


X = 1 . 0 
D09I=1,m 
D06I J=10 . 120 

6 K t I J ) = 199 728 
D07L=15, 115,10 

7 K l L J =1 0 5520 

NN= [[ Y( I 1-YmINI 1 )]/ 1 YMAX* YM IN III 11*100 .+15,5 
K INN) =143344 

PRlNTd.X* [K[M]# m=10.120J 

B FORMA T I LH F7,1V1X,6<U,15M,15A1,15A1#15A1#15A1,15A1#15A1) 
9 X=X+1 , 0 
RETURN 
END 



APPENDIX C 


MASS ATTRACTION OF THE TILT METER 

The level in the mercury pools of the tilt meter changed with the posi- 
tions of the attracting mass. Because the distances between the mercury pools 
and the attracting lead mass are large, the effect is very small. The mercury 
surface of the tilt meter pool is always vertical to the_acceleration vector con- 
sisting of the earth gravity g and the mass attraction A of the lead assembly. 
The tangent to the surface in pool number i has an angle a * to the horizontal 
plane and the tangent to the surface in pool number 2 has an angle a 2 (Fig- C-l) 
The vectors A 1 and A 2 have x- and y-components. 



A 2X a u 


Figure C-l . Mass attraction of the mercury in the tilt meter. 
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Since a i , a 2 , Ajy, and A 2 y are very small, 


and 


4x 


«l = 


*2X 


CV 2 = 


A good approximation of the imaginary equipotential surface is a sphere. The 
line connecting the two pools is the secant line of the circle which is the inter- 
section curve of the sphere with the x-y plane. The angle a between the secant 
line and the x-axis is the geometrical correlation 

q- = ~ + o ' 2 ) rad. 

This angle a is the output value of the tilt meter ; a was computed for the three 
positions of the attracting mass. 

Position up a u = 0. 680x1 0“ s rad 

Position middle o' m = 0.595xi0~ 8 rad 

Position down a ^ = 0. 496x1 0 -8 rad . 

The values also agree with those computed according to the potential theory. 
The movement of the mass causes a tilting of the mercury surface. The input 
tilting is almost a sine wave with the amplitude of a = 0. 1±0. 07x1 0 -8 rad. 

The data analysis of the tilt meter data AY(I) shows an amplitude at the fre- 
quency FO of 

BY(25) = (0.3±0.2) 10"® rad. 


The value has a large standard deviation and is in the same magnitude as the 
computed amplitude for mass attraction. Therefore, no tilting of the test pad 
connected with the movement of the attracting mass was large enough to 
influence the calibration of the MESA in the laboratory. 
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